Does the Sun Shine by pp or CNO Fusion Reactions? 
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We show that solar neutrino experiments set an upper limit of 7.8% (7.3% including the recent 
KamLAND measurements) to the fraction of energy that the Sun produces via the CNO fusion cycle, 
which is an order of magnitude improvement upon the previous limit. New experiments are required 
to detect CNO neutrinos corresponding to the 1.5% of the solar luminosity that the standard solar 
model predicts is generated by the CNO cycle. 



In 1939, Hans Bethe described in an epochal paper 0| 
two nuclear fusion mechanisms by which main sequence 
stars like the Sun could produce the energy corresponding 
to their observed luminosities. The two mechanisms have 
become known as the p — p chain and the GNO cycle . 
For both the p — p chain and the CNO cycle the basic 
energy source is the burning of four protons to form an 
alpha particle, two positrons, and two neutrinos. Thus 
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*He +2e+ + 2^^ + < 25 MeV (to the star). (1) 



In the p ~ p chain, fusion reactions among elements 
lighter than A = 8 produce a characteristic set of neu- 
trino fluxes, whose spectral energy shapes are known but 
whose fluxes must be calculated with a detailed solar 
model. In the GNO chain, with ^^C as a catalyst, ^^N 
and '^^O beta decays are the primary source of neutrinos. 

The first sentence in Bethe's paper reads: "It is shown 
that the most important source of energy in ordinary 
stars is the reactions of carbon and nitrogen with pro- 
tons." Bethe's conclusion about the dominant role of the 
CNO cycle relied upon a crude model of the Sun. Over 
the next two and a half decades, the results of increas- 
ingly more accurate laboratory measurements of nuclear 
fusion reactions and more detailed solar model calcula- 
tions led to the theoretical inference that the Sun shines 
primarily by the p — p chain rather than the GNO cycle. 
Currently, solar model calculations imply |^ that 98.5% 
of the solar luminosity is provided by the p — p chain and 
only 1.5% is provided by CNO reactions. 

In recent years, there have been many analyses of so- 
lar neutrino oscillations, essentially all of which assumed 
that the CNO neutrino fluxes were equal to their pre- 
dicted standard solar model values (see Refs. |§, |[ ||] 
and references cited therein). However, from the earli- 
est days of solar neutrino research, a primary goal of the 
field was to test the solar model prediction that the Sun 
shines by the p ~ p chain and not by the GNO cycle . 



This goal has largely been ignored in the last decade or so 
as solar neutrino experiments concentrated on the more 
accessible, higher-energy *B neutrinos. In this paper, we 
return to the question of how well we can measure, or set 
an upper limit to, the CNO neutrino fluxes. 

Unfortunately, the standard solar model prediction for 
the GNO fluxes is difficult to test. Radiochemical exper- 
iments with chlorine and gallium |||, ^ |ll| do not 
measure the energy of the neutrinos detected; they mea- 
sure the rate of neutrino induced events above a fixed 
energy threshold. The neutrino-electron scattering ex- 
periments, Kamiokande and Super-Kamiokande p3| , 
provide information about neutrinos but only those that 
have energies well above the maximum energies of the ^^N 
(E^^^ = 1.2 MeV) and ^^O {E^^^ = 1.7 MeV) neutrinos. 
The radiochemical experiments are only sensitive to elec- 
tron type neutrinos, and the neutrino-electron scattering 
experiments are primarily sensitive to electron type neu- 
trinos. The heavy water experiment, SNO ||^, |l^, mea- 
sures higher energy neutrinos. The goal of uniquely iden- 
tifying CNO neutrinos is made even more difficult by the 
fact that neutrino oscillations can change in an energy 
dependent way the probability that electron type neutri- 
nos created in the Sun reach the Earth as electron type 
neutrinos [|5[ |l6| . 

Because of these complications, it was possible to find 
neutrino oscillation solutions in which 99.95 % of the 
Sun's luminosity is supplied by the CNO cycle [pT| . 
These 'large GNO' oscillation solutions describe well all 
of the measurements from the chlorine SAGE 
GALLEX ||l^, and Kamiokande solar neutrino ex- 
periments. Modern solar models do not predict a large 
CNO contribution to the solar luminosity, but the goal 
is to test experimentally-not just assume-this prediction. 

In this paper, we use data from the chlorine, SAGE, 
GALLEX, GNO, Super-Kamiokande, and SNO so- 
lar neutrino experiments, and from the recent Kam- 
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LAND iQ reactor measurements, to set an experimental 
limit on the CNO contribution to the solar luminosity 
that is an order of magnitude more stringent than the 
previous best limit 1 17 . Although individual experiments 



do not constrain well the CNO fluxes, a global solution to 
all the available neutrino data provides a powerful upper 
limit. We also discuss how well future experiments can 
do in detecting the CNO neutrinos. 

Here is our strategy. For each value of the CNO 
luminosity fraction, Lqno/ Lq, we search the two- 
component neutrino oscillation parameter space with a 
dense mesh corresponding to the neutrino mass differ- 
ence IQ-i^eV^ < Am^ < IQ-^eV^ (721 mesh points) and 
mixing angles 0.0001 < tan^ 9 < 10 (401 mesh points), as 
well as the solar neutrino fluxes (see below). [We verify 
later that our approximation of only two neutrinos does 
not limit the validity of the upper bound we derive. See 
discussion following Eq. (||).] We calculate the global 
by fitting to all the available data. 
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We carry out a global analysis IJ] of the solar neutrino 
data letting the neutrino fluxes be free variables and us- 
ing data from 80 measurements: 44 data points from 
the Super-Kamiokande zenith-angle energy distribution, 
34 data points from the SNO day-night energy spec- 
trum 1^, 1^, and 2 radio-chemical rates from CI [@| and 



Ga |, |10|, y. We define the 3a upper limit for the CNO 
neutrino fluxes by determining when — Xmin + 9 after 
marginalizing over the oscillation parameters Arn^ and 
tan^ 9 and over the other neutrino fluxes. 

Using the data provided in Ref. we calculate the 
positron spectrum in the KamLAND detector with the 
procedures described in Refs. ||l^, |l^, In the ab- 

sence of neutrino oscillations, we find (in agreement with 
Ref. |l^) 86.8 expected neutrino events above 2.6 MeV 
visible energy for the stated experimental conditions. 
The positron energy spectrum that we calculate is in ex- 
cellent agreement with the energy spectrum presented 
by the KamLAND collaboration. Further details of our 
analysis of the KamLAND and solar data can be found 
in Ref. |2|. 

We impose the 'luminosity constraint' on the solar neu- 
trino fluxes, i.e., we require that the sum of the thermal 
energy generation rates associated with each of the solar 
neutrino fluxes be equal to the solar luminosity ||2^ . The 
fraction of the sun's luminosity that arises from CNO 
reactions can be written as [E3| : 
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(3) 



where the constant is the energy provided to the star 
by nuclear fusion reactions associated with the i*'* neu- 
trino flux, Oi is the ratio of the neutrino flux <i>i(BPOO) 
of the standard solar model to the characteristic solar 



photon flux deflned by L0/[47r(A.?7.)2(lOMeV)], and 0, 
is the ratio of the true solar neutrino flux to the neutrino 
fluxes predicted by the BPOO standard solar model 
Ref. ||2^ presents a detailed derivation of Eq. (||) and the 
numerical values for the coefficients and Oi. 

We treat as free parameters all the solar neutrino fluxes 
that are normally reported in solar neutrino calculations. 
There are then 10 free parameters: the two oscillation 
parameters, Am^ and tan^ 0,and the 8 neutrino fluxes, 
p — p, pep, ^Be,^B, and hep (from the p — p chain) and 
13n 1^0, and "F (from the CNO cycle). To speed up 
the calculations, we made some approximations that we 
have checked do not affect the accuracy of our search. 
Two of the solar neutrino fluxes, hep and ^^F, are small 
as a result of nuclear physics considerations. In the initial 
search calculations, we set hep equal to its solar model 
value and ^^F = [(i^F/^N) , ^ ,1 " N. We checked 

LV ' / soiar modelJ 

that our results are unchanged if the hep solar model flux 
is multiplied by eight (present experimental bound from 
the high energy bins at Super-Kamiokande [p4[) or if we 
set the ^^F flux equal to the ^^N flux. Also, the ratio of 
the pep neutrino flux to the p — p neutrino flux is fixed to 
high accuracy because they have the same nuclear matrix 
element. We have set the ratio equal to the standard 
solar model value and have checked that our results are 
unchanged if this ratio is varied by 10% (an enormous 
change). We set the "'^'^N flux equal to the ^^O flux, which 
is expected in the limit that the CNO contribution to 
the luminosity is dominant. We also verified that the 
results of our search are unchanged if we set the ratio of 
^^N to ^^O neutrino fluxes equal to the standard solar 
model value, the ratio expected if the p — p contribution 
is dominant. Finally, we checked several intermediate 
values of this ratio to see that the upper limit we quote 
here is robust and valid in all cases. 

We find the minimum value of for each assumed 
value of the CNO luminosity fraction by marginalizing 
over the neutrino oscillation parameters and over the 
non-CNO neutrino fiuxes. We performed the calculations 
in two stages: first using only the solar neutrino data 
and second using both the solar neutrino and the Kam- 
LAND data. We carried out calculations for oscillations 
to purely active neutrinos, to purely sterile neutrinos, 
and to active-sterile admixtures as described in Ref. |2^ 
(see also last reference in Ref. jl^). We considered sterile 
admixtures that range from the maximum allowed by the 
recent KamLAND data, sin^ 77 = 13% to 25%, 50%, 
75%, as well as the extremes of 0% and 100%. For all 
values of the CNO luminosity fraction, the minimum 
was, as expected, achieved for purely active oscillations. 

Figure H summarizes our main results. The figure 
shows Ax as a function of the CNO luminosity frac- 
tion when only solar neutrino data are used (denoted by 
dotted curves) and when solar and KamLAND data are 
used (denoted by solid curves). The minimum value of 
X^, relative to which Ax^ is measured, is reached in both 
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FIG. 1: Experimental bound on Lcno/Lq. The figure shows 
how the fit worsens as one increases the assumed fraction 
of the Sun's luminosity that arises from CNO reactions. The 
dotted lines were computed using just solar neutrino data; 
the solid lines use both solar neutrino experiments and the 
KamLAND results. The curves labeled LMA were calculated 
for the favored large mixing angle MSW solution; the curves 
labeled non-LMA were calculated for the best fit of the LOW, 
SMA, vacuum, and sterile oscillation solutions. The arrows 
indicate the predicted 1.5% CNO luminosity from the stan- 
dard solar model and the ~ 8 % (see Eq. ^) 3a upper bound 
(1 dof) allowed by the chlorine, gallium, Super-Kamiokande, 
and SNO solar neutrino data |, |, |^, ^ Q and the 
KamLAND reactor data jl^. 



cases for a zero value of the CNO flux. However, as is 
apparent from Figure the global is essentially flat 
for all values of Lcno/Lq < 5%. Current experiments 
are not sensitive to CNO neutrino fluxes that correspond 
to less than 5% of the solar luminosity. For example, the 
Ax^ is only 0.5 between the best fit at Lcno/Lq — 0.0 
and the standard solar model value of 1.5%. 

In all cases, the best fit is a LMA (large mixing angle) 
MSW solution In principle, new solutions might 

have been found by allowing Lcno/Lq to be a free pa- 
rameter. In practice, no preferred new solutions are 
found for any value of Lcno/Lq. 

For completeness, we explored the entire oscillation pa- 
rameter space specified just preceding Eq. (|^), including 
the parameter ranges of all previously recognized oscil- 
lation solutions ^, ||]. The best-fit results for these 
searches for other oscillation cases are labeled 'non-LMA' 
in Figure |l|. For Lcno/Lq less than about 5%, the 



best-fit non-LMA solution is the well-known LOW so- 
lution. For the larger values of Lcno/Lq shown in Fig- 
ure |l|, the best-fit LMA solution is a vacuum solution with 
Am^ = 7.9 X IQ-^eV^ and tan^ 9 = 0.22. We conclude, 
using both solar neutrino and KamLAND experimental 
data (using just solar neutrino data), that 
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< 7.3% (7.8%) at 3(7 . 



(4) 



The recent KamLAND measurements reduce the upper 
limit by 0.5%. 

The order of magnitude improvement between the pre- 
vious hmit of 99.95%|]l^ and the present hmit, Eq. (|), is 
due to the Super-Kamiokande and SNO measurements. 
The earlier large CNO oscillation solution was confined 
to small mixing angles, SMA, which cannot fit simultane- 
ously the flat recoil energy spectrum measured by Super- 
Kamiokande [l^ and the total event rates measured by 
Super-Kamiokande and SNO ||, |l|. 

The maximum CNO neutrino flux allowed by the ex- 
isting experiments is 



(/-CNO, max < 3.41 X 10^" Cm ^ 



{Lcno/Lq) , (5) 



where 4>CN0, max = 0(^^N)max = </'(^^0)inax- 

We have verified that the upper limits given in Eq. (^) 
and Eq. (|^) are not affected by the approximation of as- 
suming that there is just one mass scale (i. e., two neu- 
trinos). We repeated the analysis assuming the standard 
three-neutrino mixing scenario invoked to explain both 
solar and atmospheric data and assumed values for ^13 
values below the CHOOZ bound H, tan2(6ii3) = 0.0, 
0.03 and 0.06 (the CHOOZ bound at 3cr Q). The min- 
imum was, as expected, achieved for tar? (6*13) — 0.0. 

New solar neutrino experiments are required to mea- 
sure the CNO contribution to the solar luminosity. 

How much can a future ^Be neutrino-electron scatter- 



ing experiment, BOREXIN O |28[ or KamLAND [iS), im- 
prove the limit given in Eq. I^l We find an approximate 
answer to this question by computing a global x^ iH" 
eluding the existing solar neutrino data, the KamLAND 
reactor data, and a simulated BOREXINO rate mea- 
surement (simulations guided by Ref. ||2^). We assume 
that the BOREXINO rate will be consistent with the 
predicted best fit point from the solar plus KamLAND 
global fit with a total error of 10% (5%) for the rate 
measurement. If these assumptions are valid, one will 
be able to either measure Lcno/Lq or conclude that 
Lcno/Lq < 5.6% (4.9)%. 

In order to measure the CNO contribution at the 1.5% 
level predicted by the standard solar model, one must be 
able to distinguish the continuum ^^N and "'^^O neutrinos 
from the ^Be and pep neutrino lines, as well as from all 
the sources of background. The appropriate analyses of 
proposed low energy neutrino-electron scattering detec- 
tors have not yet been carried out, so one cannot say for 
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sure whether or not this will be possible. But, it seems 
very difficult. The energy resolution required to measure 
the energy of the CNO neutrinos and determine their 
flux, may, however, be within the reach of low-energy 
CC experiments p9[ . 

The solar model predictions for CNO neutrino fluxes 
are not precise because the CNO fusion reactions are not 
as well studied as the p — p reactions Q and because the 
Coulomb barrier is higher for the CNO reactions, imply- 
ing a greater sensitivity to details of the solar model. 
For the standard solar model CNO neutrino fluxes, the 
1(7 errors vary between 17% and 25% A measurement 
of the CNO neutrino fluxes would constitute a stringent 
test of the theory of stellar evolution and provide unique 
information about the solar interior. 
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